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Abstract 

We consider Vachaspati's primordial magnetic field which is generated at the 
electroweak phase transition. Assuming that either the gradients of the Higgs 
field or, alternatively, the magnetic field itself are stochastic variables with a 
normal distribution, we find that the resulting magnetic field has an rms value 
in the present-day universe which is fully consistent with what is required for 
the galactic dynamo mechanism. 



1 Introduction 



The magnetic fields observed in galaxies (of the order 10^® G) can be understood as 
an amplification by a dynamo effect [||] of a weak seed field of order 10^^*^ G on a co- 
moving scale of 100 kpc. The intergalactic plasma has a large electrical conductivity, 
and the magnetohydrodynamical equation 

9B/(9t = Vx (vxB)-Vx (r^VxB) , (1) 

where t] is the inverse electrical conductivity, then implies that the magnetic lines of 
force are essentially "frozen into" the fluid. The magnetic flux through any contour 
moving with the plasma is thus constant. The collapse of the plasma into a galaxy 
enhances the magnetic field (by a factor ~ 10^), whereas the remaining necessary 
enhancement is due to the differential rotation and turbulent motion of the plasma in 
the galaxy. 

Eq. (|l]) is homogeneous in B, so if initially one has no field it follows that a field 
can never be generated. This is basically the reason for the need of a seed field. 

It has often been speculated that the seed field is of primordial origin, which means 
that it should be explained by features relevant for the early universe. Electromag- 
netism first occurs when the standard electroweak SU{2) ® U{1)y theory is broken 
down to f/(l)em- It is therefore particularly attractive that Vachaspati has explained 
the origin of a primordial field in terms of the cosmological boundary condition that 
all physical quantities should be uncorrelated over distances greater than the horizon 
distance. Since we start with the group SU{2) ® U{1)y before the electroweak phase 
transition, the resulting electromagnetic field can be constructed in a way which is 
different from the usual F^i. = dfj^Aiy — d^A^. The result is 

where is the Higgs field. At the electroweak phase transition the correlation length 
in the broken phase is ~ l/m\Y (assuming that the Higgs mass is comparable to mw). 
The field Fij is thus constant over a distance ~ 1/mw, but it varies in a random way 
over larger distances in order to respect causality. The vector Vi is also random, of 
course. Its variation is due to the fact that the field (p makes a random walk on the 
vacuum manifold of (p. The problem then is to estimate the field Fij over a length scale 
~ N/mw li N = 1, then it follows that on dimensional grounds Fij ~ ~ 10^^ 
G, with probably an uncertainty of ±1 in the exponent 0. For large, one should 
use a statistical argument. In p| it was argued that the gradients are of order 1/a//V, 
since makes a random walk on the vacuum manifold with A0 ~ V^, and since 
Ax ~ N. Thus Vi is, in a root mean square sense, of the order and hence 

Fij is of order Taking further into account that the flux in a co-moving circular 

contour is constant, the field must decrease like l/a(t)^, where a(t) is the scale factor 
[0]. Using the fact that in the early universe a goes like the inverse temperature, the 
field was then estimated to behave like 

{F^,)T ~ ^ (3) 

when the temperature of the universe is T. For a scale of 100 kpc this leads to 
{Fij)now ~ 10"'^'^ G, which is far too small to explain the galactic fields (unless there 
exists some large scale amplification mechanism). 



In the present paper we shall present a different statistical scenario where the 
gradient vectors are taken to be the basic stochastic variables. Our considerations have 
been inffuenced by the fact that in large scale 3- dimensional computer simulations of 
the dynamo effect the computer uses a random initial magnetic field configuration. 
We find that the very interesting expression (H) obtained in can be interpreted 
statistically in such a way that the mean magnetic field satisfies 

{F,,)t = , ~ ^ • 

Comparing the root-mean-square value with Eq. (^), we observe that the scaling 
behavior is weaker by ^/N. This means that for a scale of 100 kpc 

sJ{F^j)tod.y ~ 10-''G , (5) 

which is very close (if not equal) to the value desired for the dynamo effect. 

It should be emphasized that in estimates it is reasonable, from the point of view 
of the dynamo effect, to calculate the rms value (or, more precisely, to estimate the 
rms value of the projection of the random field on the dynamo eigenfunctions) [^], 
whereas in actual computer simulations the field configuration should be generated by 
the computer subject to the two conditions (^). 



2 Discussion of the averaging procedure 

Let us now turn to the detailed arguments. We wish to consider random fields walking 
around in space in a certain number of steps. Thus we replace the continuum by a 
lattice, where the points are denoted by greek letters a, . . .. We want to estimate the 
magnetic field over a linear scale (which at most is equal to the horizon scale). Thus, 
we consider a curve consisting of N steps in the lattice, and we define the mean value 

— 1 ^ 

^ = T7E^"' > (6) 

i=i 

where B is a component of the magnetic field, and where the lattice points are on 
the curve. 

Now this curve is arbitrary, and we could take any other curve. We therefore define 
the average (...), which averages over curves spanning an A^^ lattice, i.e. over all space 
(this is well defined on a lattice space; one could e.g. take the set of all curves that 
are parallel to one of the sides of the lattice). Then, for example, 

— 1 ^ 

i=i 

which means that for each curve with steps the mean value B is computed, and this 
is done for a set of curves which span an A^^-lattice, and the average is then computed. 
Therefore {B) depends in general on N, but for simplicity of notation we shall leave 
out the explicit reference to this dependence. We wish to emphasize that the ensemble 
average (0) takes into account the field value at each lattice point, so that the average 
is really over the whole lattice volume. 

Similarly, one can define higher moments such as 



together with quantities hke 

Note that in (|]) the sum is over curves of length steps of the non-local quantity 



3 Random Higgs gradients 



In this section we shall present our main statistical assumptions. The general point 
of view is that when one has a random system it is necessary to specify the statistical 
distribution and also which variable is to be considered as a stochastic variable. These 
two specifications are the necessary boundary conditions. 

In the stochastic variable was taken to be the Higgs field itself which varies 
over the vacuum manifold. However, it is clear that also the gradient vectors Vi are 
stochastic, and in our scenario we assume that these vectors are the relevant stochastic 
variables. This is because they directly specify whether there is a magnetic field or not, 
whereas this is only true indirectly for the Higgs field itself. Also, the vectors Vi are 
relevant for questions of alignment between neighbouring domains. Thus, a scenario 
which takes the gradient vectors as the basic stochastic variables is rather natural. 
In this scenario the vectors and the resulting magnetic field have only short-range 
correlations. 

We now return to the expression (|^) of the magnetic field in terms of the Higgs 
gradients Vi. It is convenient to split these fields into real and imaginary parts, 

\^,(x) = i?,(x)+^/i(x) , (9) 

where Ri and Jj are real vectors. We consider the system at a fixed time. The 
cosmological boundary condition is then that Ri and /j are random fields. We now 
make the following assumptions: 

(i) The random fields have a Gaussian distribution. Thus, the mean value of some 
quantity Q is given by 

{Q) = l[l 'IK^Q ^-x[Rt-m?~MiT-{h)? ^ (10) 

where D is a normalization factor defined such that (1) = 1, and A is a measure 
of the inverse width. The quantities Ri and Jj are the mean values of Ri and 
li defined along a curve of length steps.Q Thus, eq. ([I0|) is relevant for a 
3-dimensional world which is an A^^ lattice. 

(ii) We assume that the mean values are isotropic, i.e. = (-R2) = (-R3) and 
(/i) = {h) = (Is). 

Assumption (i) is certainly the simplest way of implementing lack of correlation 
of the gradient vectors over distances compatible with the horizon scale, whereas as- 
sumption (ii) is natural as there is no reason to expect any preferred direction. 



It should be noted that the distribution (10) factorizes into an /-part and an im- 



part. Thus, for any expectation value consisting of J's and i?'s one has factorization, 

(^■■■K ^•••^) = (^n---^D(^---^) • (11) 

This property turns out to be very useful in computing the higher moments. 

^Note that this imphes that the mean value in a point is assumed to be equal to the mean value 
computed along all curves of length N. Thus the mean values can depend on N. 



4 The expectation value of the magnetic field 

First we consider the expectation value of a component Bi of the magnetic field. From 
the expression (||) we find that 

Bi = -^ijkFjk = —i SijkVjyk = '^^ijkRjlk ■ (12) 

Thus 

1 AT T N 

B, = j,T.BT = '^e,,,-Y.RTiT ■ (13) 

Hence 

r, AT 

i=i 

2 ^ 

= J^^MH - {RiWT - ih)) + N{R,){h)) . (14) 
1=1 

Now, due to the factorization (0), the first term on the right-hand side of the last 
Eq. ([T^ ) vanishes,^ and hence 

(Bj) = 2e,ik(Ri)(Ik) = Sjik {(Ri)(Ik) - (Rk)(Ii)) = (15) 

because of the isotropy assumption (ii). Consequently the mean value of the magnetic 
field vanishes, as announced in the Introduction. 

It should be noticed that if we did not assume isotropy, then (Bj) 7^ in general. 
If we then use Vachaspati's argument 0, and (/«) behave like I/a/ZV, and hence 
from (^) one would find (Bj) ~ 1/N. Thus, our method of averaging over curves 
leads to the same result as found previously, if we do not assume isotropy. However, 
we believe that isotropy is a natural assumption. 



5 The root mean square of the magnetic field 

The second order moment is given by 



af3 



Y^{{R^Rf'){rif') - {RtR^){I^If)} 



(16) 



Q/3 



where we used the factorization (|TT]). Now 



= n 

7 

1 

2A 



D 



+ n 



D 



{R,)R'^+{R,)R'^-{R,){R,) 



(,17) 



and similarly for (I^'Ij). Further we have e.g. 



n 



^^^^ T>Pp-KR^-{R)f 
D ' 



n 

7 



d^W 
D 



?/3 m \^„-A(i^T-(i^))2 



(R^ - {R,))e 



^Because (i?"' — {Ri)) = (/^' — {Ik)) = for symmetry reasons. 



i.e., the mean value in a given arbitrary point (3 on the lattice is equal to the mean 
value computed over all curves. Using and ([T8|) in (|^) we get 



{b]) = ^,i:{-^,+\{{i?+{W))+^,^^^^ ■ (19) 

The first term is 0{^) = 0{j^). The last term, being the square of the mean value, 
actually vanishes because of isotropy: If we take 

(R) = l={r,r,r) ; (!) = l={c,c,c) ; {W = ; (If = , (20) 

then 

{W(I)' - mm = r'c' - (^(3rc))' = . (21) 
Thus we conclude that the rms value of the magnetic field is given by 



as announced in the Introduction. 

The reason for this slow decrease is the fact that isotropy prevents the mean value 
from entering in {Bi) and {B ), and that the correlations of the gradient vectors are 
of short range. 

6 Another probability distribution 

It should be mentioned that the scenario developed in the previous section is by no 
means unique from the point of view of producing short range correlations. Instead of 
assuming that the vector Vi has a random distribution one could make the assumption 
that the magnetic field B^ itself has a random distribution, with the probability 

Y[l[e-^iK->"LL. , (23) 

i=lj = l ^ 

where D is a normalization factor. Here we have assumed that = 0. Then one 
has 

(^rsj) = (24) 

and hence 

1 ^ 1 

= Jf. UiB'-'f) - , (25) 

i.e. the same as the previous result. It should be noted that this result appears in 
spite of the fact that the distribution ( p3D is very different from the distribution (p!l|) 
of the vector field, since e.g. (^) contains correlations between the li and Ri fields. 
Also, it should be noticed that the reason for the result ( P3| ) is that the magnetic field 
has only short-range correlations. 



2 ^ / 3 1 1 



The distribution (^) is the one which is usually assumed in solid state physics 
when dealing with a random magnetic field. In the continuum version it reads 

e-^Jd'^B^i^f-^Bix) (26) 

and one then has 

(5.(x)5,(y)) = ^M'(x-y) , (27) 

where the (5-function is assumed to be smeared. 

Let us also comment on the case when (B) ^ 0. For the fiuctuations we would 
have 

(^') = 4 E(^"^") = ^? + 4 E((^" - (^))(^' - (^))) • (28) 

a/3 afi 

With only short range correlations, i.e. with 

((5- _ (5))(5'« - (5))) = 5°^((5" - (5))2) (29) 

this gives 

(fi') = (fi)^ + 0(1) . (30) 

Thus, if one has {B) ~ 0(1/A^) as in ref. 0, and if the correlations are of short range, 
then the dominant term is the fiuctuations (9(1/A^) in Eq. (|30|). Therefore, even in 
this case one has (-B)^ ^ (-B ), and hence the field should be estimated from the 
rms- value \/ {B^) when is large, not from (5), and the rms- value is again effectively 



of order l/yN . 

7 The energy-momentum tensor 

We should briefly discuss the consequences of the results described in Sections 4-6 
for the energy- momentum tensor T^^. After a calculation analogous to that in Eqs. 
(16)-(22) we obtain 



1 / —2, , ,^2, 



((/,,)(/,) + . (31) 



Contraction over i and / of course reproduces Eq. (0). We see that the expectation 



value in Eq. (plD has two terms, an isotropic and an isotropic term. This is also 



carried over to the energy-momentum tensor 



^00 = 2"^^' "^^^ ^ 2^*'"^^ ^ BiBi . (32) 
Consider flrst the case (//) = {Ri) = 0. Then 

{Tu) = ^Su{B') - {B,Bi) = ^uiB') , (33) 

so that the equation of state is still isotropic: p = ^p. However, in the case where (/;) 
and (Ri) are no n- vanishing we obtain instead of Eq. ( ^2]) 



mi) = j^Su + ^{{U){Ii) + {R^){Ri)) . (34) 



Hence in this case {Tu) has an additional anisotropic part. 

If the field Bi itself is a Gaussian random field, we find that 

(Too) = \{{B - {B)f) + \{Bf , 

mi) = \6u[l{{B-m + {W]-mm , (35) 

so in general there also appears an anisotropic part in {Tn). If (B) is 0{1/N) as in 
the anisotropic part is subdominant because of Eqs. ( pBD - (0). Thus, to the leading 
order {Tn) is isotropic. 

In principle the magnetic energy-momentum tensor discussed above contributes to 
the expansion of the universe. However, in the present universe its effect is extremely 
small and can safely be ignored. In the early universe this may not be so. Our re- 
sults show that in principle there may exist non-isotropic stresses of magnetic origin, 
which in a statistical sense could contribute to produce turbulence in the primeval 
plasma, giving possibly rise to an early "universal" dynamo effect. Magnetically gen- 
erated plasma flows might also be important e.g. for the dynamics of the QCD phase 
transition. 



8 Consequences of the magnetic field 

We assume that at the electroweak scale T = To ~ 100 GeV the coherence length 
of the rms fleld is .^o — l/^o so that in terms of the physical distance L, we have 
N = L/^Q. The magnetic fleld is frozen at that time, so that at later times the original 
coherence length is redshifted by the expansion according to 

m = ^^o . (36) 
ao 

The frozen-in magnetic field is also redshifted by the expansion of the universe. Thus 
at later times at the distance scale L, 




B,Ut,L) = Bo ^ ^ = 5o ? t P^ , (37) 



where Tgto = 0.301 Mp/ ^J'gjTo) with g^: the effective number of degrees of freedom, 
and U ^ 1.4 x 10^(^0/^^)"^ yrs is the time when the universe becomes matter dom- 
inated; for definiteness, we shall adopt the the value floh'^ = 0.4, which is the upper 
limit allowed by the age of the universe. We shall also assume that Bq ~ 10^"^ G. 

We may easily find from Eq. (|37| ) the size of the cosmological field today, which 
could have acted as the seed field for the dynamo mechanism. Taking t = 1.5 x lO^'' 
yrs and L = 100 kpc (corresponding to = 1.0 x 10^^), we find that today the 
cosmological field at the scale of intergalactic distances is 

Brms = 4 X 10"^^ G . (38) 

This seems to be exactly what is required for the numerical dynamo simulations to 
produce the observed galactic magnetic fields of the order 10~^ G. The inherent un- 
certainties in the estimate ( pSj) are: the value of Qoh"^ used for computing t*; the time 
at which the magnetic field froze, or Tq; the actual value of the field Bq. Therefore 
one should view (153) as an order-of-magnitude estimate only. 



We should also check what other possible cosmological consequences the existence 
of the random magnetic field, Eq. (|37|) , may have. Let us first note that the energy 
density ps in the rms field is very small. In the radiation dominated era we find that 
the energy density within a horizon volume V is 

1 r^'H 3 / T \ ^ 1 

The horizon distance is th = 2t so that ps ~ T^/Mp <^ p^, and the magnetic field 
contribution to the total energy density is negligible. 

In principle, magnetic fields could modify primordial nucleosynthesis. For instance, 
it has been argued that protons actually become heavier than neutrons in a large 
enough magnetic field. In our case, however, the magnetic field is glued to the charges 
according to Eq. (|l|) so that the relative velocity is zero. Note that this prevents the 
charged particles in the plasma from accelerating by emitting synchrotron radiation, 
as would happen if there were a constant background field. At the onset of nucleosyn- 
thesis, at about T ^ 1 MeV, the effect of the rms field on the weak reaction rates 
which change protons to neutrons and determine the crucial n/p-ratio turns also out 
to be negligible. Since at T ~ 1 MeV n ^ p reactions have scattering lengths of the 
order of the horizon length, from Eq. (^) we find that the relevant field at that scale 
is only Srms — 1500 G. Also, creating a thermal population of right-handed neutrinos, 
disastrous for the successful prediction of primordial light element abundances, via 
scattering of left-handed neutrinos off the magnetic field may not be possible in 
our case because the mean squared length of the field fiuctuation is expected to be very 
short, 0{1/T). This means that the ul i^r transition probability would be very 
much suppressed 0. However, this issue can only be settled by detailed dynamical 
considerations, for example by computer simulations. 

Perhaps more interesting is the role of the rms field at the QCD phase transition 
[|]. It is believed to be of first order, but the details of bubble nucleation depend on 
the largely unknown dynamical details of QCD. In simplistic nucleation theory one 
obtains, by comparing the nucleation rate with the Hubble rate, for the size of the 
critical bubble about 10 fm. The distance d between nucleation centers depends on 
the amount of supercooling, but with reasonable assumptions (i>10~^ m. This is then 
the scale at which the bubbles of new phase will feel the random background field. We 
find 

5rms(tQCD,rf)<1.8x 10^ G . (40) 

Whether this has an effect on the QCD phase transition or not depends to a large 
extent on whether the quarks are glued to the flux lines. We shall not discuss this 
issue further. 

9 Conclusions 

We have shown that if the derivative of the (logarithm of the) Higgs field is a random 
variable, then Vachaspati's construction leads to a magnetic field in the present day 
universe which has the right order of magnitude from the point of view of the galactic 
dynamo mechanism 0]. 

It should also be mentioned that this result is more general, as is clear from Sect. 6. 
This is because if by some mechanism a Gaussian random magnetic field is generated 
at the electroweak phase transition (where it is always of order over a correlation 
length on dimensional grounds), then Eq. (^5]) shows that the rms- value behaves like 



l/\/iV, if {B) = 0. However, even if (5) 7^ 0, it follows from Eq. (^) that the rms 
value is at least of order I/a/ZV. This is because either (-B)^ is larger than or equal to 
1/iV, in which case the rms value is larger than or equal to or (-B)^ is less than 

in which case the fluctuation term of order dominates, and the rms- value 
is of order 1/\/N. Thus, if there exists any mechanism which produces a Gaussian 
random magnetic field at the electroweak phase transition, it will produce a result 
which is larger than or equal to the field required by the dynamo effect. It therefore 
appears that there exists a good case for the primordial origin of the observed galactic 
magnetic fields. 

If we had considered a magnetic field performing a random walk in a given vol- 
ume, or a random magnetic flux through a given surface, we would have obtained 
the scalings and 1/iV, respectively. Then one still would have had to weight 

these spatial averages by the statistical distribution, but this would then have induced 
double counting of the lattice points. Therefore we have considered random walk only 
along a given curve, together with an ensemble average over all curves spanning a 
given volume. Therefore our method takes into account the field value at each lattice 
point, and in this sense the ensemble average is really over the whole lattice volume. 
Only a detailed dynamical simulation can however tell for sure whether our seed field 
actually gives rise to the observed galactic magnetic field. A definite prediction of our 
scenario is that the seed field has {B) = by virtue of the averaging procedure, and 
that only the rms field is non-vanishing. 
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